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Satellite trace gas sensing modeling 
for missions such as ASCENDS 


• Simulation framework to predict performance of trace gas sensing in the atmosphere from space 

• Spectroscopic (excitation wavelengths optimizations) and lidar analysis 

• Applicable to active sensing missions such as ASCENDS and a variety of molecules including C0 2 , CH 4 , N 2 0 etc. 

• Projected open source distribution 


ASCENDS mission overview 
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Identify anthropogenic 
generation of C02 sources 
and sinks to enable 
effective carbon trading 



Requirements and approach 

• ASCENDS is identified as a medium size Tier II NASA mission in the NRC Decadal Survey and is currently slated for 
2019 launch 

• ASCENDS will deliver laser based remote sensing measurements of Global C0 2 mixing ratios 

• (XC0 2 ) to a precision of 0.3 percent on horizontal scales of 100-km over land and 200-km over oceans and passive 
CO measurements for C0 2 interpretation. 

• ASCENDS is the logical extension of OCO and GOSAT 


Benefits for climate science 

• Quantify global spatial distribution of atmospheric C0 2 on scales of weather models 

• Quantify global spatial distribution of terrestrial and oceanic sources and sinks of C0 2 during day/night over all 
seasons. 

• Provide a scientific basis for future projections of C0 2 sources and sinks through data-driven enhancements of 
Earth-system process modeling. 


NASA LaRC ASCENDS approach 

• < 0.3% accuracy (~ lppm) in C0 2 mixing ratio resolution is required NASA Langley Research Center (LaRC) is 
developing an intensity modulated continuous wave (IMCW) laser absorption spectrometer based remote sensing 
scheme for the detection of C0 2 at 1.57 microns and 0 2 at 1.26 microns from space based platforms 

• Multiple wavelengths with differential absorption (DIAL) technique are used 1.26 micron band for 0 2 sensing is 
selected to obtain surface pressure 

• 1.26 micron band provides architectural and spectroscopic advantages 

• For our experiments two candidate wavelengths in this band will be utilized. Nominally, wavelengths around 1.262 
and 1.271 microns have been identified in initial tests lines in these two sub-bands are being further analyzed for 
sensitivity to environmental parameters 

Modeling framework for integrated 
path space lidar performance 

estimates 

Fig l.l Previously under development 








Multi-wavelength LIDAR 
simulation program 

(Space-based and airborne 
LIDAR sensing) 


Shell GUI program to 

run the AER LBLRTM 


Line-by-line transmission 
calculation program 

1) HITRAN with Voigt profile 

2) Custom lineshape routines 


HITRAN database 
management program 

(Analog of JAVA HAWKS) 
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Fig 1.2 New additions 


/ 

5 Additional program 
to estimate temperature 
and water vapor 
interferences 

(Based on pre-analyzed 
satellite data) 
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Temperature sensitivity 

(As function of wavelength) 
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Water vapor interferences 

(As function of wavelength) 
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Alternative methods 

(As function of wavelength 
combinations) 
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Load data from ft* 


Pre-analyzed atmospheric data for 

error analysis 

Analysis and processing of the MERRA dataset 
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Fig 2.1 Pre-analyzed MERRA 
atmospheric data 


Analysis of pressure, temperature, 
and humidity fields averaged for 
all Earth coordinates for a span of 
year 2009 (the frequency of each 
combination is color coded on a 
logarithmic scale of 100 to 10 7 ). 

Grid: 

1 Kelvin x 1 mbarx 0.0001 kg/kg 
Temperature x Pressure x Humidity 


Fig 2.2 Altitude distribution 
□f the analyzed MERRA 
fields 

Analysis of MERRA data at 
different pressure levels 
Showing density of temperature 
and humidity combinations at 
selected altitudes 


Specific humidity, kg/kg 


Annual MERRA data for all locations is analyzed to establish the density of points within each 
temperature x Pressure x Humidity grid bin 

Combining individual layer error estimates 


Fig 2.3 Averaged temperature 
and pressure profiles 

MERRA atmospheric parameters 
of pressure and temperature as 
a function of layer number for the 
analyzed fields on model layers 
(inst6_3d_ana_Nv) dataset;. 

(a) - Center positions of the total 
pressure levels as a function of layer 
number, (b) Annual average of the 
temperature values for year pressure 
level as a function of layer number. 
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Fig 2.4 Individual layer and 
cumbinad total path 
analysis approach 

Overview of separate layers 
analysis methodology based 
on pre-analyzed MERRA 
atmospheric fields. 



Global annual temperature sensitivity 
analysis for CD 7 and 0 7 bands 


Global annual temperature sensitivity for alternative C0 2 bands 
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Fig 3.1 TamparaturE 
sensitivity analysis for the 
1.57 micron C0 2 band 

Total absorption spectrum (a) 
and the temperature induced 
error (b) for the 1.57 pm band 
of C0 2 due to annual variations 
in atmospheric temperature for 
year 2009. 


Wavenumber, cm ' 1 



Fig 3.2 Temperature 
sensitivity analysis for the 
2.05 micron C0 2 band 

Total absorption spectrum (a) 
and the temperature induced 
error (b) for the 2.05 pm band 
of C0 2 due to annual variations 
in atmospheric temperature for 
year 2009. 


Wavenumber, cur 1 

Preliminary annual cumulative error due to annual global temperature variations 

Global annual temperature sensitivity for alternative 0 2 bands 


Fig 3.3 TemperaturE 
sensitivity analysis for 
the 0 2 A-band 

Total absorption spectrum (a) 
and the temperature induced 
error (b) for the A-Band of 
0 2 due to annual variations 
in atmospheric temperature 
for year 2009. 
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Fig 3.4 TempEraturE 
sensitivity analysis for 
The I.2G-I.27 micron 0 
band 


Total absorption spectrum (a) 
and the temperature induced 
error (b) for the 1.26-1.27 pm 
band of 0 2 due to annual 
variations in atmospheric 
temperature for year 2009. 
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Analyzed MERRA grid data is used in line-by-line calculations for each layer to get coefficients of variation In 
each layer. Such individual layer coefficients of variation are scaled with layer molecular density and added 
together to get the total path relative uncertainty. 


Temperature sensitivity is comparable for the A-Band and 1.26-1.27 micron 0 2 bands. The 1.26-1.27 micron 
band allows achieving lower temperature sensitivities due to closer proximity to line centers. 
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Wavelength instability effects 

Fig 4.1 Laser wavelength jitter analysis example for selected C0 2 spectral lines 
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Fig 4.2 Laser wavelength jitter analysis example for selected 0 2 spectral lines 








Laser wavelength stability requirements in different spectral regions of C0 2 and 0 2 molecules to achieve an accuracy of 0.21% 
in relative optical depth accuracy 


Molecule 

Center 
investigating 
the band 

Ref. 

Possible lines and wavelengths 

Approximate 
Maximum 
optical depth 

Maximum total jitter for < 0.21% 
optical depth accuracy 

Wavenumber, 

cm' 1 

Wavelength, 

mn 

MHz 

pm 

cm' 1 

0 2 

GSFC 

1,2 

13075.819 

764.7705 

0.87 

16 

0.031 

0.00053 

0 2 

GSFC 

1,2 

13075.944 

764.7631 

1.47 

11 

0.022 

0.00037 

0 2 

GSFC 

1,2 

13076.064 

764.7561 

3.15 

8 

0.016 

0.00027 

0 2 

GSFC 

1,2 

13077.297 

764.684 

0.46 

1400 

2.73 

0.047 

0 2 

LaRC 

2 

7841.9866 line 

1275.1871 

1.95 peak 

1 

0.0054 

0.000033 

0 2 

LaRC 

2 

7920.596 

1262.5312 

~ 0.68 

130 

0.69 

0.0043 

co 2 

GSFC 

2 

6359.9674 line 

1572.335 

~1.9 peak 

1 

0.0082 

0.000033 

co 2 

LaRC 

2 

6364.922 line 

1571.111 

- 1.38 peak 

1 

0.0082 

0.000033 

co 2 

JPL 

additional 

example 

4875.6905 

2050.992 

2.6 

2 

0.028 

0.000067 

co 2 

JPL 

3 

4875.6265 

2051.018 

0.9 

4 

0.056 

0.00013 

co 2 

JPL 

4 

4875:882 

2050.911 

0.89 

4 

0.056 

0.00013 


Notes: Wavelength positions are approximate based on the data previously reported. The values are presented for specific frequencies or entire 
lines (marked as "line"). Intensities are specified for either exact wavelength potions or peaks of the lines. The laser jitter is specified for the total 
wavelength variation span of the laser wavelength. 
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Conclusions and further work 

• Framework applicable to any molecule 

• Framework operated in the Beta version, further debugging and improvements are ongoing 

• Simulation programs are used for the analysis of temperature sensitivity, wavelength stability effects and the 
weighting functions calculations 

• Further analysis and reduction of water vapor interferences is under investigation 

• Alternative measurement approaches are being studied 

• Verification of the modeling framework components is planned 

• The advantages of the implemented methodology is the reduction in required number of calculations (ordinary PC 
may be used) 
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